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The nucleon-nucleon interaction is constructed by means of the J-matrix version of inverse scattering theory. 
Ambiguities of the interaction are eliminated by postulating tridiagonal and quasi-tridiagonal forms of the potential 
matrix in the oscillator basis in uncoupled and coupled waves, respectively. The obtained interaction is very accurate 
in reproducing the NN scattering data and deuteron properties. The interaction is used in the no-core shell model 
calculations of 3 H and 4 He nuclei. The resulting binding energies of 3 H and 4 He are very close to experimental values. 



I. INTRODUCTION 

The nucleon-nucleon (NN) interaction is one of the 
most important ingredients in microscopic nuclear struc- 
ture studies. Ideally, the NN interaction should be de- 
rived from the quark-gluon theory. The modern status of 
QCD, however, makes it possible to predict reaction cross 
sections only at high enough energies; the QCD-based 
derivation of the NN potential describing the nucleon- 
nucleon interaction at the low energies that are of pri- 
mary importance for nuclear physics applications is im- 
possible at the moment. 

Nucleon-nucleon potentials conventionally referred to 
as 'realistic', are derived from the meson exchange the- 
ory. Modern realistic NN potentials like Bonn Ar- 
gonne Q, Nijmegen ,3|, etc., are carefully fitted to the ex- 
isting experimental data on NN scattering and deuteron 
properties. Unfortunately, none of the known NN inter- 
actions provides a completely satisfactory description of 
the trinucleon and other light nuclei. To overcome this 
deficiency, meson exchange or phenomenological Q 
three-nucleon forces are usually introduced. Impressive 
progress was achieved recently in the description of the 
trinucleon and 4 He binding energies with realistic NN 
and three-nucleon forces ((j- However, the three-nucleon 
force parameters in such studies are sometimes fitted to 
the trinucleon binding and some of them may not be con- 
sistent with the parameters of the two-body interaction. 
In one very detailed study, when the three-body interac- 
tion parameters were chosen consistently with the two- 
body parameters, the three-nucleon force contribution to 
the triton binding energy was shown to be negligible • 
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Impressive progress using effective field theory has re- 
cently been reported (see review in Ref. The versions 
that provide the most accurate fit to the nucleon-nucleon 
properties Q use a momentum-space cutoff and are still 
quite strong at short distances. The match with the nu- 
clear many body model space cutoff is unclear; additional 
renormalization is required for typical model spaces that 
are feasible. We aim in this paper to have high qual- 
ity descriptions of the phase shifts with softer potentials 
whose cutoff is well-matched to the anticipated applica- 
tion in many-body systems. 

Various microscopic models have been designed for the 
studies of few-body systems. It was demonstrated in 
Ref. [T(i| that all modern realistic microscopic models pro- 
vide approximately the same results for the 4 He ground 
state. The no-core shell model [n|, is one of these 
models. This model can be used not only for the few- 
body nuclear applications but also, with modern com- 
puter facilities, for microscopic studies of heavier nuclei 
with the number of nucleons A up to A ~ 12 12]. The 
no-core shell model is based on a wave function expansion 
in a many-body oscillator function series with the aim to 
describe bound states and narrow resonances treated as 
bound states. 

The oscillator basis matrices of the modern realistic 
NN potentials are very large and cannot be directly used 
without a severe truncation in the many-body no-core 
shell model calculations. As a result, the convergence 
of the calculations appears to be slow. This deficiency 
is conventionally addressed by constructing the so-called 
effective NN interaction (see, e. g., [ll|). Ideally the 
effective NN interaction should reproduce in the finite 
model space the results of the infinite model space cal- 
culation. In a realistic application, the construction of 
the effective NN interaction is a complicated problem 
involving various approximations. 

In this contribution, we construct the NN interaction 
by means of the J-matrix version of inverse scattering 
theory [H Q IH- The matrix of the NN potential 
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in the oscillator basis is obtained for each partial wave 
independently. Therefore, in our approach we derive the 
NN interaction as a set of potential matrices for different 
partial waves. We reproduce the experimental NN scat- 
tering data and deuteron properties with small potential 
matrices. Our NN interaction can be treated as an effec- 
tive interaction since its matrix can be directly used in 
the no-core shell model calculations without additional 
truncation. However, our effective NN interaction re- 
produces the energy spectrum and other observables in 
a many-body system as well as deuteron properties and 
NN scattering data. From this point of view, our NN 
interaction can be treated as a realistic one as well. Our 
interaction is not related to the meson exchange theory, 
however we shall see that we obtain the deuteron and 
scattering wave functions that are very close to the ones 
obtained with realistic meson exchange potentials. 

The potential derived by the J-matrix inverse scat- 
tering approach is ambiguous. The ambiguity originates 
from the phase-equivalent transformation suggested in 
Ref. (see also 0, 113 and references therein) . The 
ambiguity is eliminated in the present approach by a phe- 
nomenological ansatz that the potential matrix in the 
uncoupled partial waves is tridiagonal. Therefore our 
potentials are Inverse Scattering Tridiagonal Potentials 
(1STP). The non-central nature of the NN interaction 
is manifested in the coupling of some partial waves, and 
the tridiagonal potential ansatz should be extended to 
allow for the coupling of these partial waves. We pos- 
tulate phcnomenologically the simplest generalization of 
the tridiagonal form of the potential matrix in this case; 
however, we refer to our potentials as 1STP in the cases 
of both uncoupled and coupled partial waves (though, 
strictly speaking, it is not correct in the later case), ft is 
just the tridiagonal ansatz that brings us to the scattering 
wave functions which are very close to the ones provided 
by the meson exchange realistic NN potentials. How- 
ever, in the case of the coupled sd waves we perform a 
phase equivalent potential transformation to improve the 
description of the deuteron properties. 

The suggested ISTP are used in the no-core shell model 
calculations of 3 H and 4 He. We shall see that the pre- 
dicted 3 H and 4 He binding energies are very close to the 
experimental values. We do not use three-nucleon inter- 
actions, yet our predictions of the 3 H and 4 He bindings 
are approximately of the same accuracy as the predic- 
tions based on the best realistic meson exchange two- 
nucleon plus three-nucleon forces. 

Here we would like to mention some recent papers 
where other approaches to the problem of constructing 
high-quality effective interaction were utilized. The au- 
thors of Ref s. [lfi l2fT| added phenomenological non-local 
terms to a cut-off Yukawa tail of the realistic NN poten- 
tials. The obtained interaction reproduces the 3 H bind- 
ing energy. The additional non-local terms do not reduce 
the rank of the potential energy matrix in the oscillator 
basis of the underlying realistic NN interaction. There- 
fore the use of this interaction in the shell model studies 



requires the construction of the shell model effective in- 
teraction. 

A very interesting approach is the construction of the 
low momentum NN potential Vi ow -k from the realistic 
NN interactions (see the review in Ref. H3]). The use of 
Vio W -k in the shell model applications still requires the 
construction of the shell model effective interaction but 
this problem is simplified. The effective interaction ob- 
tained from Vi ow -k was used successfully in various shell 
model applications (see, e. g., H^)- Unfortunately it is 
still unclear whether this interaction provides the correct 
binding of three-body and four body nuclear systems. 
Contrary to Vi ow -k, our ISTP is designed for the direct 
use in the shell model applications. 

The paper is organized as follows. In the next Section 
we present the single channel J-matrix inverse scattering 
approach, derive ISTP in the uncoupled partial waves 
and discuss their properties. The derivation and discus- 
sion of the ISTP properties in the coupled partial waves 
can be found in the Section III. The results of the 3 H and 
4 He calculations are presented in the Section IV. A short 
summary of the results can be found in Section V. 



II. SINGLE CHANNEL J-MATRIX INVERSE 
SCATTERING APPROACH AND ISTP IN 
UNCOUPLED NN PARTIAL WAVES 

The J-matrix formalism in the quantum scatte ring 
theory was initially proposed in atomic physics |23| . 
Within the J-matrix formalism, the continuum spec- 
trum wave function is expanded in an infinite series of 
L 2 functions. This approach was shown to be one of 
the most efficient and precise methods in calculations of 
photoionization 0, HH HH an d electron scattering by 
atoms |27). In nuclear physics the same approach has 
been developed independently [2^, |2i| as the method of 
the harmonic oscillator representation of scattering the- 
ory. This method has been successfully used in various 
nuclear applications allowing for the two-body contin- 
uum, e. g. nucleus- nucleus scattering has been studied 
in the algebraic version of RGM based on the J-matrix 
formalism (see the review papers H^El); the effect of 
A and neutron decay channels in hypernuclei production 
reactions has been investigated in Refs. [33, H3, etc. The 
approach was extended to the case of true few-body scat- 
tering in |34| and utilized in the studies of the monopole 
excitations of the 12 C nucleus in the 3a cluster model 
in Ref. [3^. It was also used in the studies of double- A 
hypernuclei in Ref. |36| and of weakly bound nuclei in the 
three-body cluster model in Refs. plill7lll8|. 

The J-matrix version of the inverse scattering the- 
ory was suggested in Refs. [lj, [lj, Ula]. The discus- 
sion of the general formalism below follows the ideas of 
Refs. mm in, however, some formulas are presented 
here in a manner that should be more convenient for the 
current application. The tridiagonalization of the inter- 
action obtained by the inverse scattering methods have 
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not previously been discussed in the literature, hence the 
corresponding theory and results are new. 

The oscillator-basis J-matrix formalism is discussed in 
detail elsewhere (see, e. g., [23, We present here 

only some relations needed for understanding the inverse 
scattering J-matrix approach. 

The Schrodingcr equation in the partial wave with or- 
bital angular momentum I reads 



H l ^ lm (E,r) = E^ lm (E,r) 



The wave function is given by 



1 



*i m (E,r) = - Ul (E,r)Y lm (r), 



(1) 



(2) 



where Yi m (r) is the spherical function. Within the J- 
matrix formalism, the radial wave function ui(E,r) is 
expanded in an oscillator function series 



ui 



(E,r)=^a nl (E) R nl (r), 



(3) 



71=0 



where 



R nl (r) = (-1)' 



2n! 



r r(n + I + 3/2) \r 



i+i 



r 2 \ / r 2 



L"(x) is the associated Laguerre polynomial, the oscil- 
lator radius ro — ^Jh/mui, and m is the reduced mass. 
All energies are given in the units of the oscillator basis 
parameter hui. 

The wave function in the oscillator representation 
cini(E) is a solution of the infinite set of algebraic equa- 
tions 



(H l nn ,-6 m ,E)a nn (E) = 0, 



(5) 



where the Hamiltonian matrix elements H l nn , = T l nn , 
V* t , the kinetic energy matrix elements 
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rtn-i = ~^n{n + l + 1/2), 



T^ n = -(2n + l + 3/2), 



(6a) 
(6b) 



Tl n+1 = _Iv(n + l)(n + / + 3/2), (6c) 

and the potential energy V within the J-matrix formal- 
ism is approximated by the truncated matrix with ele- 
ments 



y-l - J Kin' ir n an d n ' < N; 
nn ' ~ ' if n or n' > N. 



(7) 



In the inverse scattering J-matrix approach, the potential 
energy is constructed in the form of the finite matrix of 
the type (J2J); therefore the J-matrix solutions with such 
an interaction are exact. 

In the external part of the model space spanned by the 
functions Q) with n > N, Eq. j£j) takes the form of a 
three-term recurrence relation 



+ (T l nn -E)a nl (E) 



^n, n+1 On+l, 



i(E)=0. (8) 



Any solution of Eq. (JHJ) is a superposition of the fun- 
damental regular S n i(E) and irregular C n i{E) solu- 
tions |23ll37l|. 

a nl (E) = cos 6(E) S nl (E) + wo. 5(E) C nl {E), (9) 



where 

Sm(E) = 



it ro n 



T(n + l + 3/2) 



(10) 



C nl (E) = (-1) 1 



it ro ni 



T(n + l + 3/2) r(-/ + l/2) 



x exp (-1.J _ I - 1/2, -/ + 1/2; q 2 ), (11) 

<l>(a, b; z) is a confluent hypergeometric function |38| . q = 
\/2E, and 5(E) is the scattering phase shift. 

The wave function in the oscillator representation 
o-ni(E) in the internal part of the model space spanned by 
the functions (0J with n < N, can be expressed through 
the external solution on + i,i(E): 



a n i(E) = ^ n N T N N+1 ON+l,l(E). 
The matrix elements, 

N 



E 

A'=0 



Ey-E 



(12) 



(13) 



are expressed through the eigenvalues E\ and eigenvec- 
tors (n\X) of the truncated Hamiltonian matrix, i. e. E\ 
and (n\X) are obtained by solving the algebraic problem 



£ H l nn ,(n'\\) = E x {n\X), 



n < N. 



(14) 



The matrix element Sf/vjv is of primary importance in the 
calculation of the phase shift 5(E): 



t&n 5(E) 



Sni (E) — &nn T N ^ N+1 Sn+i, 1 (E) 
' Cm (E) - % NN T l N ~ C N+1J (E) • 



(15) 



In the direct J-matrix approach, we first solve Eq. 1141) 
and next calculate the phase shift 5(E) by means of 
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Eq. JT3J. In the inverse scattering J-matrix approach, 
the phase shift 5(E) is taken to be known at any energy 
E and, instead of solving lfl4")l. we extract the eigenvalues 
E\ and the eigenvectors (n\X) from this information. 

First we assign some value to N, the rank of the de- 
sired potential matrix [see Eq. 0]. Generally, with a fi- 
nite rank potential matrix it is possible to reproduce the 
phase shift 5(E) only in a finite energy interval; larger 
A supports a larger energy interval. However, from the 
point of view of many-body applications, it is desirable 
to have A as small as possible. 

The components a n i (E) of the wave function in the os- 
cillator representation, should be finite at arbitrary en- 
ergy E. This is seen from Eqs. (fT2*)> - l(T3j) to be possible 
at the energies E = E\, A = 0, 1, ... , A only if 



a>N+ij(E\) — 0. 



(16) 



Knowing the phase shift, we can calculate cln+ij(E) at 
any energy E using Eq. 0. Therefore we can solve nu- 
merically the transcendental equation (|16(l and find the 
eigenvalues E\, X — 0, 1, ... , A. 
Due to Eq. (Q3J), 



o.n+i,i(E) 



\E-E X ), 



where 



da N+ ij(E) 



dE 



(17) 



(18) 



E—E\ 



Now it is easy to derive from Eqs. H12|) (|13|l the following 
equation: 



a Nl (E x ) = \(N\X)\ 2 a? T l N<N+1 , 



(19) 



hence we should have 



A' 



]T(A|A)(A|A) = 1. 



(22) 



A=0 



or, equivalently, 



Generally the set of (A| A) obtained by means of Eq. I|2U|) 
violates the completeness relation <|22[) . Therefore this 
set of (-/V|A) ideally describing the phase shifts, cannot be 
treated as the set of last components of the normalized 
eigenvectors (n|A) of any truncated Hermitian Hamilto- 
nian matrix; in other words, the set of (A|A) violating 
Eq. ill' 21) cannot be used to construct a Hermitian Hamil- 
tonian matrix. 

To overcome this difficulty, we fit Eq. (|22|l by changing 
the value of the component (A|A = A) corresponding to 
the highest eigenvalue E x= n. This modification spoils 
the description of the phase shifts 5(E) at energies E 
different from E\, A = 0, 1, ... , N. We restore the 
phase shift description in the energy interval [0, E x= n-i] 
by variation of E\=n ■ From the above consideration it is 
clear that larger A values make it possible to reproduce 
phase shifts in larger energy intervals [0, E\=n—i]- 

There is an ambiguity in determining the potential 
matrix describing the given phase shifts 5(E): any 
of the phase equivalent transformations discussed in 
Refs. [la d Ell [see also Eqs. below] that 

do not change the truncated Hamiltonian eigenvalues E\ 
and respective eigenvector components (A|A), results in 
a potential matrix that brings us to the same phase shifts 
5(E) at any energy E. Additional model assumptions are 
needed to resolve this ambiguity. As was already men- 
tioned, we assume the tridiagonal form of the potential 
matrix. We now discuss the construction of the tridiago- 
nal potential matrix supposing A and the sets of E\ and 
(A|A) to be known. 

If the potential matrix is tridiagonal, the equa- 
tions H14|) can be rewritten as 



IW>| 2 



ClNl(E\) 



(20) 



Within the J-matrix formalism, both a^i(E) and 
a N+i.i(E) fit Eq. and can be calculated using this 
equation at any energy E. Hence, one can also calcu- 
late a£ by means of Eq. (|18|l . Therefore the components 
(A|A) can be obtained from Eq. (|20H (the sign of the 
components (A|A) is of no importance). 

Equations l|16l) and (|20() provide the general solution 
of the J-matrix inverse scattering problem: solving these 
equations we obtain the sets of E\ and (A|A), and these 
quantities completely determine the phase shifts 5(E). 
However (A|A) are supposed to be the components of the 
eigenvectors (n\X) of the truncated Hamiltonian matrix 
[see Eq. (|14|) ] that should fit the completeness relation 



i4(0|A}+<(l|A)=£ A (0|A>, 



H 1 n .n-i(N - 1|A) + H l NN (N\X) = E X (N\X). 



(23a) 



Kn-i(n - 1|A) + H l nn (n\X) + < n+1 <n + 1|A) 

= E x (n\X) (n = l, 2, ... , N-l), (23b) 



(23c) 



The unknown quantities in Eq. 



are the component 



(A — 1 1 A) and the Hamiltonian matrix elements H l N N _ 1 
and H l NN . We multiply Eq. l(23c|l by (X\N), sum the 
result over A, and use the completeness relation Q21H to 



obtain the formula for the calculation of H l NN : 



11 NN 



N 

E 

A=0 



E X (N\X) 2 



(24) 



JV 



5>|A)<A|n') 



(21) 



A=0 



The Hermitian conjugate of Eq. (|23c|) reads: 

(A | AT - l)H l N N ^ + (X\N)H l NN = (X\N)E X . (233c) 
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We multiply Eq. I|23c[l by Eq. I|23cl-c[l , sum the result over 
A, and use the completeness relation <|21|) to obtain the 
following expression for the calculation of H l N N _ 1 : 



H 



AT, N—l 



\ 



JY 



J2El(N\Xr-(H 



NN) 



(25) 



A=0 



Generally, the sign in the right-hand-side of Eq. 125|) is 
arbitrary. Here we use an additional assumption that 
the off-diagonal Hamiltonian matrix elements H l n n±1 are 
dominated by the kinetic energy so that the sign of these 
matrix elements is the same as the kinetic energy matrix 
elements T^ n±1 [see Eqs. ©]■ This assumption brings 
us to the minus sign in the right-hand-side of Eq. I|25(l . 

Now equation H23c[l can be used to calculate the last 
unknown quantity, 



(JV - 1|A) 



f 



TTl 

JV, N — l 



(E X (N\X) - H l NN (N\X)) . (26) 



We now turn to equation (|23b[) with n = N—l. 
This equation contains one more term than Eq. 123cfl . 
however this term does not include unknown quantities. 
We perform with Eq. I|23b|) exactly the same manipu- 
lations to obtain expressions for H l N _ 1 N _ 1; H l N _ 2 JV _ 1 
and (N-2\X). Setting n = N-2 in Eq. (|23rj|l . we obtain 
the expressions for H l N _ 2 N _ 2 , H l N _ 3 N _ 2 and (N — 3|A), 
etc. Equation H23afl is needed only to calculate the last 



matrix element Hq . As a result, we obtain the following 
generalization of Eq. (|2"4"|) valid at n = N, N — 1, ... , 0: 



TABLE I: Non-zero matrix elements in huj units of the 8huj 
ISTP matrix in the j sq partial wave. 



V 1 

1 el r, 



V 1 = 

v ra, ra + 1 *ra + l, n 



-0.3706925910512869 

1 -0.1599160886224698 

2 0.1395932055925835 

3 0.2668242073073204 

4 0.0414909332158313 



0.1340546812405571 
0.01647436916961609 
-0.1334461921366397 
-0.07869019612934602 



TABLE II: Non-zero matrix elements in hu units of the Thuj 
ISTP matrix in the pi partial wave. 



V, 



V, 



I 

ra, ra + 1 



V, 



I 

ra + 1, n 



0.10619936477245400 

1 0.32183202739863150 

2 0.38227890301930240 

3 0.08818666274780007 



-0.09441150969281098 
-0.19861423056402480 
-0.12529300192170380 



TABLE III: Non-zero matrix elements in huj units of the 8ftw 
ISTP matrix in the 1 d,2 partial wave. 



v; 



v, 



I 

n, ra + 1 



V, 



I 

ra+1, n 



-0.04182464628865646 

1 -0.1129604626451339 

2 -0.1276115098155470 

3 -0.02554669840505408 



0.03831247883572108 
0.06873518464832284 
0.04042212068310880 



H 



x 

E 

A=0 



E x {n\Xf 



(27) 



The equations 



H 



n, n — l 



N 



j2Ei{n\xr-(HLy-( H ^ n+1 y 

\ A=0 



and 

(n-l|A) 



(28) 



1 / N 



-^ m <«|A)-< Il+1 (n + l|A)J (29) 

are valid at n = N — 1, N — 2,... , 1. Equations lj23)) - 
(|29|l make it possible to calculate all unknown quantities. 
After calculating the Hamiltonian matrix elements H l nn , , 
we derive the ISTP matrix elements by the obvious equa- 
tions 



V 1 = H 



I 

nn 



Vn, nil — Hn, nil 



± n, nil- 



(30a) 
(30b) 



The above theory is used to construct the NN ISTP 
matrix elements in uncoupled partial waves. We use as 
input the np scattering phase shifts reconstructed from 
the experimental data by the Nijmegen group 3]. The 
oscillator basis parameter hoj — 40 MeV. Usually in the 
shell model calculations, the complete xtuv model space 
is used, i. e. all many-body oscillator basis states (con- 
figurations) with tti < k where the single-particle 
state oscillator quanta = 2ni + U, are included in the 
calculation. Thus, to be applicable to all p-shell nuclei 
in accessible model spaces, we suggest the 8ftiu and 7huj 
ISTP, i. e. the rank of the ISTP matrix N is chosen so 
that 2N + I = 8 in the partial waves with even orbital 
angular momentum I and 2N + 1 = 7 in the partial waves 
with odd orbital angular momentum I. 

The non-zero matrix elements of the obtained ISTP in 
uncoupled partial waves are presented in Tables IjT lyTTTl 
(in hu = 40 MeV units). 

In Figs. ITlJTol we present the results of the phase shift 
and scattering wave function calculations with our ISTP 
in the uncoupled partial waves. The phase shifts are seen 
to be better reproduced by ISTP up to the laboratory en- 
ergy i?i a b = 350 MeV than by one of the best realistic 
meson exchange potentials Nijmegen-II. Some discrepan- 
cies are seen only at large energies. These discrepancies 
can be eliminated by using larger N values. This is il- 
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TABLE IV: Non-zero matrix elements in hui units of the 7 hui 
ISTP matrix in the 1 f$ partial wave. 



n VL 

' " r Tin 


y! _ yl 

'n,n + l *n + l,ri 


0.04238710037363047 


-0.02790556099208952 


1 0.07474001110555983 


-0.02815383549650696 


2 0.02511618088977574 




TABLE V: Non-zero matrix elements 


in hui units of the 7hui 


ISTP matrix in the ^po partial wave. 




n V 1 


V 1 , , — V 1 , , 


—0 13674752057396140 


01511502604716686 


1 0.08786870226069166 


-0.10590497118041760 


2 0.23624887864971810 


-0.08040102075340183 


3 0.04909915603358606 




TABLE VI: Non-zero matrix elements 


in hui units of the 7hu> 


ISTP matrix in the pi partial wave. 




n V n7l 


V 1 — V 1 


0.08893328127606703 


-0.09288011075116410 


1 0.33899943058663640 


-0.21111518227376910 


2 0.36158649481733020 


-0.09828565221995666 


3 0.05167268571103811 




TABLE VII: Non-zero matrix elements in hui units of the 8hu> 


ISTP matrix in the 3 d,2 partial wave. 




n V nn 


Vn, n + 1 — n 


-0.20024057805517500 


0.1193321938724192 


1 -0.28898789873267020 


0.1463047726434774 


2 -0.25522202901437920 


0.0792277802117809 


3 -0.05421394437761595 




TABLE VIII: Non-zero matrix elements in hui units of the 


7hu> ISTP matrix in the 3 fs partial wave. 


n VL 


Vn, n + 1 — + n 


0.02629214811765302 


-0.01394097030190910 


1 0.03463672270715756 


-0.01259217885072639 


2 0.01119624135204766 





lustrated in phase shifts of odd partial waves presented 
in Figs. E| El ED ED These are the results of the 
phase shift calculations with the 9hu> ISTP in addition 
to the Ihui ISTP phase shifts. It is interesting that the 
differences between the 7Hu> ISTP and 9hu> ISTP wave 
functions in odd partial waves are too small to be seen 
in Figs. 0] El El El Eleven at large energies. We note 
also that the use of 7hui ISTP instead of 9Huj ISTP in 
the 3 H and 4 He calculations, result in negligible differ- 
ences of the binding energies, wave functions, etc. The 
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FIG. 1: 1 so np scattering phase shifts. Filled circles exper- 
imental data of Ref . Q jsolid line — realistic meson exchange 
Nijmegen-II potential Q phase shifts; dashed line — ISTP 
phase shifts. 
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FIG. 2: 1 so np scattering wave functions at the laboratory 
energies E lah = 2, 10, 50, 150, 250 MeV. Solid line — realis- 
tic meson exchange Nijmegen-II potential 3] wave functions; 
dashed line — ISTP wave functions. 

ISTP np scattering wave functions at different energies 
are very close to the Nijmegen-II wave functions both in 
odd and even partial waves. In other words, these ISTP 
wave functions can be regarded as realistic. 



III. TWO-CHANNEL J-MATRIX INVERSE 
SCATTERING APPROACH AND ISTP IN 
COUPLED AW PARTIAL WAVES 

In the case of the nucleon-nucleon scattering, the spins 
of two nucleons can couple to the total spin S = (sin- 
glet spin state) or to the total spin S* = 1 (triplet spin 
state). In the case of the singlet spin state, we have only 
uncoupled partial waves in the nucleon-nucleon scatter- 
ing. In the case of the triplet spin state, the total angular 
momentum j = 1 + 1 can be obtained by the coupling of 
the total spin S = 1 with the orbital angular momentum 
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FIG. 3: 1 pi np scattering phase shifts. Filled circles — ex- 
perimental data of Ref. Q; solid line — realistic meson ex- 
change Nijmegen-II potential J| phase shifts; dashed line — 
7hu> ISTP phase shifts; dotted line — 9hco ISTP phase shifts. 
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FIG. 4: 1 p\ np scattering wave functions at the laboratory 
energies E lah = 2, 10, 50, 150, 250 MeV. Solid line — realis- 
tic meson exchange Nijmegen-II potential Q wave functions; 
dashed line — 7hu> ISTP wave functions; dotted line — 9hu> 
ISTP wave functions. 
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FIG. 5: 1 d,2 np scattering phase shifts. See Fig. 0for details. 
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FIG. 6: 1 d,2 np scattering wave functions at the laboratory 
energies £ lab = 2, 10, 50, 150, 250 MeV. See Fig.glfor details. 



The wave function in the coupled waves case is 

* = £~u r (£,r)|r), (31) 



I. On the other hand, the higher triplet-spin partial wave 
of the same parity with the orbital angular momentum 
I' = I + 2, can have the same total angular momentum 
j = 1 + 1 = 1' — 1. Such partial waves are coupled due 
to the non-central nature of the NN interaction. The sd 
coupled partial waves (the coupling of the 3 s 1 and 3 d\ 
partial waves) and pf coupled partial waves (the cou- 
pling of the 3 p2 and 3 / 2 partial waves) are of special 
interest for applications. The case of the sd coupled par- 
tial waves is of primary importance due to the existence 
of the only np bound state (the deuteron). The coupled 
equations describing the NN system in the coupled par- 
tial waves, are of the same structure with the coupled 
equations describing the two-channel system. In other 
words, the description of the coupled waves in the NN 
scattering is formally equivalent with the description of 
the two-channel scattering. 



where |T) is the spin- angle wave function which includes 
the spin variables of two nucleons coupled to the total 
spin 5 = 1, the spherical function Yi rm (r), and the cou- 
pling of the channel orbital momentum lr with the to- 
tal spin 5 into the total angular momentum j; u r (E,r) 
is the radial wave function in the given formal channel 
r = {^r,j}- Generally there are two independent so- 
lutions for each radial wave function u r (E, r). To dis- 
tinguish these solutions it is convenient to employ the 
isf-matrix formalism associated with the standing wave 
asymptotics of the wave function: 

u r{ri) (E,r) — > £ (srrjJ^-Krr^E) nj^)) . 

(32) 

Here the index Tj distinguishes independent radial func- 
tions w r (r 4 ) {E : r) in the channel T, Kpp^E) is the K- 
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FIG. 7: 1 /3 np scattering phase shifts. See Fig.[3]for details. FIG. 9: 3 po np scattering phase shifts. See Fig. [3]for details. 
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FIG. 8: 1 /3 np scattering wave functions at the laboratory 
energies £i ab = 2, 10, 50, 150, 250 MeV. See Fig.Hfor details. 



FIG. 10: 3 po np scattering wave functions at the laboratory 
energies £ lab = 2, 10, 50, 150, 250 MeV. See Fig.^Jfor details. 



matrix, and j[ (x) and n ; (x) are spherical Bessel and Neu- 
mann functions. The advantage of the if-matrix formal- 
ism is that the radial functions u r ^ r .\ (E, r) defined ac- 
cording to their standing wave asymptotics i\.>'2\) are real 
contrary to the more conventional S-matrix formalism 
with complex radial wave functions which are asymptot- 
ically a superposition of ingoing and outgoing spherical 
waves. The iif-matrix K-pTi(E), of course, can be ex- 
pressed through the S'-matrix. However it is not the S'- 
matrix but the so-called phase shifts St and Sr t in each 
of the coupled partial waves T and I\ and the mixing 
parameter e that are usually published as functions of 
the energy E in the experimental and theoretical inves- 
tigations. The S-matrix can be parametrized in terms 
of Sr, <5r ; and e. However for the present application 
it is more convenient to express the if-matrix elements 
directly through <5r, Sr t and e (see Refs. [3flL liflj): 



K sd (E) = K ds (E) 



^ tan S s + tan 2 e ■ tan 8d 

1 — tan e ■ tan S s ■ tan Sd 

. . tan Sd + tan 2 e ■ tan S s 

K dd{E) = " J 1 x 1 T> 

1 — tan e ■ tan 6 S ■ tan dd 



(33a) 
(33b) 



tane 



cos S s ■ cos Sd ■ (1 — tan e ■ tan S s ■ tan S d ) 

(33c) 

To be specific, we have specified the case of the coupled 
sd waves where the channel indexes T and I\ take the 
values s or d. In the case of the coupled pf waves, one 
substitutes the indexes s and d by the indexes p and / 
in the above expressions and in other formulas in this 
section. 

Within the inverse scattering J-matrix approach, the 
potential in the coupled partial waves is fitted with the 
form: 



V =l^ l nF > V nn> (n'T'\. (34) 

r.T' "=0 n'=0 



Here V^, = (nT\V\n'V) is the potential energy matrix 
element in the oscillator basis 

\nT)=R nlr (r)\T), (35) 
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FIG. 11: 3 pi np scattering phase shifts. See Fig.[3]for details. FIG. 13: 3 d,2 np scattering phase shifts. See Fig.^for details. 
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FIG. 12: 3 pi np scattering wave functions at the laboratory FIG. 14: z di np scattering wave functions at the laboratory 
energies E Xzh = 2, 10, 50, 150, 250 MeV. See Fig.Hfor details. energies £ lab = 2, 10, 50, 150, 250 MeV. See Fig.Hfor details. 



where the radial oscillator function R n i r (r) is given by 
Eq. (J3J and |T) is the spin- angle function. Different trun- 
cation boundaries can be used in different partial 
waves r. 

The multi-channel J-matrix formalism is well known 
(see, e. g., [23l H3) and we will not discuss it here in 
detail. The formalism provides exact solutions for the 
continuum spectrum wave functions in the case when the 
finite-rank potential V of the type (f 341) is employed. In 
the case of the discrete spectrum states, the exact solu- 
tions are obtained by the calculation of the corresponding 
S'-matrix poles as is discussed in Refs. ^ n P ar ~ 

ticular, the deuteron ground state energy Ed should be 
associated with the S'-matrix pole and its wave function 
is calculated by means of the J-matrix formalism applied 
to the negative energy E = Ed- 

Within the J-matrix formalism, the radial wave func- 
tion iir(r 4 ) (E, r) is expanded in the oscillator function 
series 

oo 

a nr(ri)(E) Rni T (r). (36) 

n=0 



In the external part of the model space spanned by the 
functions i|35|) with n > Nr , the oscillator representation 
wave function a^rfroO^) n ^ s the three-term recurrence 
relation JSJ. Its solutions corresponding to the asymp- 
totics (1321) are 

«„r(ro (E) = 6rr t S nlr (E) + K rTz (E) C nlr {E). (37) 

Equation (|37|l can be used for the calculation of 
a nr(r )(E) with n > Ny if the coupled wave phase shifts 
8r and <5rv and the mixing parameter e are known. 

The oscillator representation wave function a n Y(Vi)(E) 
in the internal part of the model space spanned by 
the functions ll33|) with n < Nr, can be expressed 
through the external oscillator representation wave func- 
tions aN r +i,r(r t )(E) as 

a nr{ r z ){E) = &nN r , T Nr,,N r ,+i «jv r ,+i,r'(r,)( £ )- 

(38) 
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FIG. 15: 3 /3 np scattering phase shifts. See Fig.[3]for details. 
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FIG. 16: 3 /3 np scattering wave functions at the laboratory 
energies £ lab = 2, 10, 50, 150, 250 MeV. See Fig.Hfor details. 



The matrix elements, 



'nn' 



N 

£ 

A'=0 



{nv\\') (Ayr) 

Ey-E 



(39) 



where N = N-p + iVjv + 1, are expressed within the di- 
rect J-matrix formalism through the eigenvalues E\ and 
eigenvectors (nr|A) of the truncated Hamiltonian matrix, 
i. e. E\ and (nr|A) are obtained by solving the algebraic 
problem 

2V r , 

J2T, H nn'(n'r'\X)=E x (nT\\), n<N r . (40) 

T' n'=0 

Here H^n> = (nT\H\n'T') are the Hamiltonian matrix 
elements. 

Within the inverse J-matrix approach, we start with 
assigning some values to the potential truncation bound- 
aries Nr [see Eq. (1341) ] in each of the partial waves T. As 



a next step, we calculate the sets of eigenvalues Ex and 
respective eigenvector components (./VrT|A). This can be 
done using the set of the J-matrix matching conditions 
which are obtained from Eq. 138[l supposing n = Nr- In 
more detail, these matching conditions are (to be spe- 
cific, we again take the case of the coupled sd waves so 
the channel indexes T and I\ take the values s or d) 



a N sS (s)(E) - £ Gsr> Tx r , <Nr ,+i a Nr , +1 x'(s)(E), 
r'=s,d 

(41a) 

a N d d(s) ( E ) = £ Qdv T£ r , )A r r ,+i a>N r ,+i,T'(s){E), 

T'=s,d 

(41b) 

a N B s(d)( E ) = £ G s r> 4 r ,,jv r ,+i aN r ,+i,r'(d)(E), 

(41c) 



r<=s,, 



and 



l N d d(d) 



(E) = 



£ 

T'=s,d 



GdV n 



,+x a N r ,+i,r>(d)(E), 

(41d) 



where we introduced the shortened notation 



Srr< = S& r 



N 

£ 

A'=0 



(N r T\X') (X'lNr'V) 
Ey-E 



(42) 



To calculate a N r , r ^.(E) and a N +1 r , r JE) entering 
Eqs. I|41|l . we can use Eq. i|37|l with the if-matrix 
elements expressed through the experimental data by 
Eqs. (jSnj- Therefore Gss, Gsd, Gds and Gdd are the only 
unknown quantities in Eqs. (|41|l and they can be ob- 
tained as the solutions of the algebraic problem 1411) at 
any positive energy E. 

These solutions may be expressed as 



Gss = 



Gdd = 



Ass(E) 
A dd (E) 

T N d ,N i+ l A ( E Y 



(43a) 
(43b) 



and 



Gds 



f2EK sd 



^^N s ,N,+l ^N d ,N d + l A(E) 



(43c) 



where 
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A SS (E) = (s NsS (E) + K SS {E) C NsS {E)) (s NM (E) + K dd {E) C NM (E)^j - K 2 sd (E) C NsS {E) C NM (E), (44a) 
A dd (E) = (s Ns+hs (E) + K SS (E) C Ns+1 , s (E)) (s Ndd {E) + K dd (E) C Ndd (E)^j - K 2 sd {E) C Na+x , s {E) C Ndd {E), (44b) 

and 

A(E) = (S N3+1 , S (E) + K SS {E) C Ns+1 , s (E)) (s Nd+1 , d (E) + K dd (E) C Nd+1 , d (E)) - K 2 sd {E) C Ns+1 , s {E) C w , d (E). 

(44c) 

To derive Eq. (|43c|) , we used the following expression for the Casoratian determinant [3, : 

JC(C, S) = C n+1>l (E) S nl (E) - S n+ld (E) C nl (E) = r ° . (45) 

^ J n,Ti+l 



It is obvious from Eqs. H42[) and (|43|) that the eigenval- 
ues E\ can be found by solving the following equation: 

A(E X ) = 0. (46) 

The eigenvector components (iVrT|A) can be obtained 
from Eqs. (|43ajl - (|43b|) in the limit E — ► E\ in the same 
manner as Eq. 1|2(J[I in the single-channel case: 

l(A^|A)| 2 = Ass(J ^. A (47) 

1 N 3 ,N a + l ^ 

and 

\(N d d\X)\ 2 = (48) 

1 N d M d + l ^ 

where 

Equations H47[) - l|48|) make it possible to calculate the ab- 
solute values of (N s s\X) and (N d d\X) only. However the 
relative sign of these eigenvector components is impor- 
tant. This relative sign can be established using the re- 
lation 

(Ex) 

(N d d\X) T& d>Nd+1 ~ a Na+1As) (E x ) 

_ aN d +i, d (d)(E\) ^ 

a N a + l,s(d)(E\) 

that can be easily obtained from Eqs. H41[) . 

Using Eqs. lf4l)|) - (|5U|l we obtain all eigenvalues E\ > 
and corresponding eigenvector components (iVrr|A). For 
example, in the case of the coupled pf waves when the 
NN system does not have a bound state, all eigenval- 
ues E\ are positive and by means of Eqs. (|46[l - i|5U[l we 



obtain a complete set of eigenvalues Ex — 0, 1, ... , 
N and the complete set of the eigenvector's last com- 
ponents (iVrr|A) providing the best description of the 
'experimental' (obtained by means of phase shift analy- 
sis) phase shifts 5\{E) and 63(E) and mixing parameter 
e. However, as in the case of the uncoupled waves, we 
should take care of fitting the completeness relation for 
the eigenvectors (nT\X) that in the coupled wave case 
takes the form 

JV 

^(nr|A)(A|n'r') = <W£rr'. (51) 

A=0 

Due to Eq. (|5ip. in the two-channel case, we should 
perform variation of the components (7Vrr|A) associated 
with the two largest eigenenergies Ex=n and Ex=n-i to 
fit three relations 

JV 

^<AT ri r 1 |A)<A|iV ri r 1 )=l ! (52a) 

A=0 
JV 

^<7V riri |A)<A|iVr 2 r 2 ) =0, (52b) 

A=0 

and 

JV 

^(jV r2 r 2 |A)(A|iVr 2 r 2 ) = 1. (52c) 

This immediately spoils the description of the scatter- 
ing data that can be restored by the additional varia- 
tion of the eigenenergies Ex=n and Ex=n-\- As a re- 
sult, in the case of the coupled pf waves, we perform a 
standard fit to the data by minimizing \ 2 by the vari- 
ation of (N pP \X = N), (N p p\X = N —1), (N f f\X = N), 
(Nff\X = N — 1), Ex=n and Ex=n-i- These 6 parame- 
ters should fit three relations l|52|l , hence we face a simple 
problem of a three-parameter fit. 
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In the case of the coupled sd waves, the np system has 
a bound state (the deuteron) at the energy E d (E d < 0) 
and one of the eigenvalues E\ is negative: Eq < 0. We 
should extend the above theory to the case of a system 
with bound states. For the coupled sd waves case when 
the np system has only one bound state, we need three 
additional equations to calculate Eq and the components 
(N s s\\ = 0) and (N d d\\ = 0). 

The deuteron energy E d should be associated with the 
S-matrix pole. As it was already noted, the technique 
of the S-matrix pole calculation within the J-matrix 
formalism is discussed together with some applications 
in Refs. 0, 0|. In the case of the finite-rank poten- 
tials of the type (|34|l . one can obtain the exact value of 



the bound state energy E d and the exact bound state 
wave function by the S-matrix pole calculation within 
the J-matrix formalism. To calculate the S'-matrix, we 
use the standard outgoing-ingoing spherical wave asymp- 
totics and the respective expression for the J-matrix os- 
cillator space wave function in the external part of the 
model space discussed, e. g., in Refs. [TTt IT81I341 137| in- 
stead of the standing wave asymptotics (|32|l and respec- 
tively modified expression i|37[l for the J-matrix oscilla- 
tor space wave function. Using the expressions for the 
multi-channel S'-matrix within the J-matrix formalism 
presented in Refs. 0, EM Hi) j it is easy to obtain the 
following expressions | 14| for the two-channel S-matrix 
elements: 



Ssd ^N„N a +l ^N d ,N d + l ^N a +l,s( E ) ^1^+1,(2 C^)} ' (53a) 



Sdd - { (pN?a( E ) ~ ^ss T N a ,N s +l C n\i,s( E )) {pN d \{ E ) ~ ^dd T N d ,N d + l C< h]+l,d( E j) 



-Gsd T N„N a +l T N d ,N d +l C N^+l,s( E ) C N d \l,d( E )} ' ( 53b ) 



and 



where 



S s d — Sds 



ir a V2Eg s 
D(E) 



(53c) 



D(E) = (C$ S (E) - Gss n M1 cW +li ,(S)) - Qu T^ w ^(B)) 

™ Ssd T N„N S +1 T N d ,N d + l C N s +l,s( E ) C< N d \lA E ) ( 54 ) 
I 



and 

C^(E)=C nl (E)±iS nl (E). (55) 

We need to calculate effi (E) at negative energy E = E d 
which can be done using Eqs. (|55(l . (|1C)|> and l|ll(l where 
imaginary values of q = qd = iy // 2\E d \ are employed. 
Extension of these expressions to the complex q plane is 
discussed in Ref. [34|. 

Since we associate the deuteron energy E d with the 
S-matrix pole, from Eqs. (|53|l we have 

D(E d ) = 0. (56) 

Assigning the experimental deuteron ground state energy 
to E d in Eq. (|56|l and substituting D(E d ) in this formula 



by its expression l|54|) . we obtain one of the equations 
needed to calculate E , (N s s\X — 0) and (N d d\X = 0). 

Two other equations utilize information about the 
asymptotic normalization constants of the deuteron 
bound state and s/ d . If the S-matrix is treated as 
a function of the complex momentum q, then its residue 
can be expressed through si s and s/ d [4lL |42| : 

i Res S lrlr , = r e^Or+'r')^ ^ . (57) 
q=tqd 

(the factor ro in the right-hand-side originates from the 
use of the dimensionless momentum q). jz£ and n = — - 
are determined experimentally. Therefore it is useful to 
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FIG. 17: Structure of the ISTP matrix in the coupled pf 
waves and of the Version ISTP in the coupled sd waves. 
The location of non-zero matrix is schematically illustrated 
by solid lines. 



rewrite equations (|57ll as 

i lim (q - iq d )S ss = r ^, 2 (58a) 

and 

i lim (q - iq d )S sd = -r i]£/^. (58b) 

Substituting S ss and S sd by its expressions (|53|l — (|54|) . we 
obtain two additional equations for the calculation of Eq, 
(N s s\X = 0) and (N d d\\ = 0). 

Clearly, in the case of coupled sd waves, we should 
also fit the completeness relation l|51|). We employ the 
following method of calculation of the sets of the eigen- 
values E\ and the components (iV s s|A) and (N d d\X). The 
E\ values with A = 1, 2, ... , N — 2 are obtained by 
solving Eq. (|46|l while the respective eigenvector's last 
components (7V s s|A) and (N d d\X) are calculated using 
Eqs. (|47I) - (|5U|I . Next we perform a x 2 fit to the scattering 
data of the parameters Eq, E x= n-i, E\ = n, (N s s\\ = 0), 



(N s s\X = N-1), (N d d\X = N), (N d d\X = 0), (N d d\X = 
N - 1), and ( N d d\X = N). These 9 paramet ers fit 6 re- 
lations (JESJ, lpE|l . (|5lc|) . (jSH3l and (|58bjl . i. e. 
we should perform a three-parameter fit as in the case of 
coupled pf waves. 

Now we turn to the calculation of the remaining eigen- 
vector components (nr|A) with n < Nr and the Hamil- 
tonian matrix elements H^, with n < Nr and n! < Nr' 
entering Eq. (|40|l . The coupled waves Hamiltonian ma- 
trix obtained by the general J-matrix inverse scatter- 
ing method is ambiguous; the ambiguity originates from 
the multi-channel generalization of the phase equivalent 
transformation mentioned in the single channel case. As 
in the single channel case, we eliminate the ambiguity by 
adopting a particular form of the potential energy ma- 
trix. 

As in the case of uncoupled partial waves, we con- 
struct %hw ISTP in the coupled sd waves. Therefore 
2N V + l r = 8, or 2N S + = 8 and 2N d + 2 = 8; hence 
N s = N d + 1. In the coupled pf waves, we construct 
7Hu> and 9ftio ISTP; clearly we again have N p = Nf + 1. 
Thus the potential matrix V„£> nas the following struc- 
ture: the submatrices V^, coupling the oscillator com- 
ponents of the same partial wave are quadratic [e. g., 
(N p + 1) x (N p + 1) submatrix V^, in the 3 p 2 wave] while 
the submatrices VF^, with r ^ P coupling the oscillator 
components of different partial waves are (Nr + 1) x Nr 
or Nr x (A^r + 1) matrices [e. g., (N p + 1) x (N p ) submatrix 
V^l, coupling the 3 p 2 and 3 f 2 waves]. Our assumptions 
are: we adopt (i) the tridiagonal form of the quadratic 
submatrices V^, and (ii) the simplest two-diagonal form 
of the non-quadratic submatrices V^, with r ^ P cou- 
pling the oscillator components of different partial waves. 
The structure of the ISTP matrices in coupled partial 
waves is illustrated by Fig. El 

Due to these assumptions, the algebraic problem 14U|) 
takes the following form: 



H° s (0s\X) + H^(ls\X) + H° d (0d\X) = E x (0s\X), (59a) 

H^(0s\X) + Hft{U\\) + H dd (0d\X) + H dd (ld\X) = E x (0d\X), (59b) 

(n = l, 2, ...,AT S -1), (59c) 

H d n s n (ns\X) + H% n+1 {n + l,s\X)+ H dd n _ x (n - 1, d\\) + H%(nd\\) + H d n d n+1 (n + 1, d\X) = E x {i d\X) 

(n=l,2, ...,N d -l), (59d) 

H s n s s ,n s -i{N s - 1, s|A) + H s N s sNs (N s s\X) + H s N d Nd (N d d\X) = E x (N s s\X), (59e) 
and 

H&,n.-i{N. - Ml A) + Hf? dNs {N.a\\)+Hfa Nd _ 1 {N d - l,d\X) + H dd dNd (N d d\X) = E x (N d d\X). (59f) 
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Even though this set of equations is more complicated than the set (|23fl discussed in the uncoupled waves case, it can 
be solved in the same manner. 

Multiplying Eqs. ifS ^ -lpf ) by (N s s\X) and (N d d\X), summing the results over A and using the completeness 
relation l|51|) we obtain 



N 



H 



N 3 N 3 



H 



dd 

N d N d 



A=0 
JV 

Y,Ex(N d d\X) 2 



(60a) 
(60b) 



and 



A=0 



N 



H s N d sNd = J2Ex(N s s\X)(X\N d d). 



A=0 



(60c) 



Now we multiply each of the equations H59e[l - lj59f|) by its Hermitian conjugate and one of these equations by the 
Hermitian conjugate of the other, sum the results over A and use (|5T)l to obtain 



H 



N S ,N 3 -1 



N 



j2ei (n sS \\) 2 (Hf? sN y (Hf* Nd y, 

\ \=o 



H 



N d ,N 3 -l 



H 



N 3 ,N 3 -1 



N 



El (N s s\X) (X\N d d) - H*J sNd (h s n * bNs 



A=0 



rrdd 
U N d N d 



(61a) 
(61b) 



and 



H 



dd 

N d ,N d -l 



N 



J2 El (N d d\X) 2 (Hff dNd Y - (Hff sN f - (Hff d 



N 3 -l) 



(61c) 



As in the case of uncoupled waves, we take the off-diagonal matrix elements N ±1 and Hff N d ±i ^° ^ e dominated 
by the respective kinetic energy matrix elements N ±1 and N ±1 and therefore choose the minus sign in the 
right-hand-sides of Eqs. (|61a|) and (|61c(l . 

By means of Eqs. (|60() - (|61|l we obtain all matrix elements HYE, entering Eqs. (|59e() - (j59f|l . Using this information, 
the eigenvector components (N s — 1, s\X) and (N d — 1, d|A) can be extracted directly from Eqs. (|59c|l - ll59fll : 



(N s -l,s\X) 



H 



and 



(N d -l,d\X) = 



N 3 ,N e -l 



rrdd 

n N d ,N d -l 



E x (N s s\X) - H s N s aNg (N s s\X) - H s N d sNd (N d d\X) 



E X (N d d\X) - H d N d dNd (N d d\X) - H d N s dN3 (N sS \X) - Hff^M - l,s\X) 



(62a) 



(62b) 



Now we can perform the same manipulations with Eqs. (|59a|l - lj59dp . We take n — N s — 1, N s — 2, ... , 1 in Eq. (|59cl) 
and n = N d — 1, N d — 2, ... , 1 in Eq. (|59d|) . Equa tions l|59c|l - l)59dll are a bit more complicated than Eqs. I|59e|l - (|59f|l . 
however the additional terms in Eqs. H59cfl - i|59d|) include only the quantities calculated on the previous step. As a 
result, we obtain the following relations for the calculation of the matrix elements ■ 



N 



N 

H dd n = Y, E ^ nd \ x y 



(63a) 
(63b) 



A=0 
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and 



N 



ff£ = X>A<n|A)(A|nd>. 



A=0 



h: 



issd 
n n,n-l 



and 



H 



n.n— 1 





N 

Y,Et(nd\\f-(Ht) 2 


- (^nfn+l) 2 














1 


" AT 






_A=0 


Tjsd I IJSS 
nn \ nn 


' n nn J n n,n+l rL n,n+l 






N 










Y,El(n S \\) 2 -(H^) 2 


~ { H n,n+l) 


~ (-^nfn-l) — (Hm) ■ 




A=0 









(63c) 



Equation H63a|l is valid for n = JV S , N s — 1, ... , while equations l|63bjl - l|63c[> are valid for n — Na, Nd — 1, ... , 0. 
For the matrix elements H^ r n _ 1 we obtain 



(64a) 
(64b) 



(64c) 



Equation (|64a() is valid for n = Nd — 1, Nd — 2, ... , 1; equation (|64b|) is valid for n — Nd, Nd — 1, ... , 1, and 
equation l|64c|l is valid for n = N 3 — 1, N s — 2, ... , 1. 

The eigenvector components (n — 1, s|A) with n = N s — 1, N s — 2, ... , 1 and (n— 1, e?|A) with rt = iV<2 — 1, Nd — 2, ... , 1 
can be calculated using the following expressions: 

(n - l,s|A) = — ^(s A (n S |A) - H£{na\\) - H°» n+1 {n + 1, a\X) - ^^(n - 1, d|A) - ff£<nd|A)) (65a) 



and 



(n - 1, d\X) = -i— (£ A (nd|A) - fl* (na|A) - ff*Ux<n + 1, s\X) - H%{nd\\) - H d n d n+1 {n + 1, d\X)) . (65b) 

n,n— 1 



Having calculated the Hamiltonian matrix elements 
ff[J , we obtain the potential energy matrix elements 
by subtracting the kinetic energy. 

We recall here that we arbitrarily assigned the values 
s and d to the channel index T but the above theory can 
be applied to any pair of coupled partial waves. The only 
equations specific for the sd coupled partial waves case 
are Eqs. ijSTfl) that are needed to account for the ex- 
perimental information about the bound state which is 
present in the np system in the sd coupled partial waves. 
In equations JSt), |@TJ|-|J50J| and one can sub- 

stitute s and d by p and /, respectively, and use them 
for constructing the ISTP in the coupled pf waves. 

We construct ISTP in the coupled NN partial waves 
using as input the np scattering phase shifts and mixing 
parameters reconstructed from the experimental data by 
the Nijmegen group We start the discussion from the 
ISTP in the coupled pf waves. 

The non-zero potential energy matrix elements of the 
obtained 7Hu> pf -ISTP are given in Table IIXI (in fiui = 
40 MeV units). The description of the phase shifts S p and 



5f and of the mixing parameter e is shown in Figs. 1181 - 
I2U1 The phcnomcnological data are seen to be well re- 
produced by the 7Hu> ISTP up to the laboratory energy 
Eub ~ 270 MeV; at higher energies there are discrepan- 
cies between the ISTP predictions and the experimental 
data that are most pronounced in the 3 p2 partial wave 
(note the very different scales in Fig. ^| and Figs. 1191 - 
12011 . These discrepancies are seen to be eliminated by 
constructing the 9hui p/-ISTP. 

Generally, for the coupled pf waves, we have 4 ra- 
dial wave function components u p /p\(E,r), u p (f\(E,r), 
Uf( p }(E,r) and Uf/fl(E,r) defined according to their 
standing wave asymptotics (|32() . We present in Figs. I21h 
[3U]the plots of these components at the laboratory ener- 
gies £i a b = 2, 10, 50, 150 and 250 MeV obtained with 
the Ihw and 9hw ISTP in comparison with the respective 
Nijmegen-II wave function components. 

It is seen from the figures that the 9hw ISTP and 
Nijmegen-II 'large' (diagonal) wave function components 
u p ( p j(E,r) and Uttf\(E,r) are indistinguishable. The 
same Ihio ISTP components differ a little from those 
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FIG. 18: 3 p2 np scattering phase shifts 8 P (coupled pf waves). 
Filled circles — experimental data of Ref. [3j; solid line — re- 
alistic meson exchange Nijmegen-II potential Q phase shifts; 
dashed line — 7hu> ISTP phase shifts; dotted line — 9hu> 
ISTP phase shifts. 
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FIG. 19: 3 /2 np scattering phase shifts 5/ (coupled pf waves) 
See Fig. H^l for details. 
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FIG. 21: Large components u p / p ^(E,r) and Uftf\(E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy i5i ab — 2 MeV. See Fig. ^| for details. 



0.002 



-0.002 



3 p 2 - 3 / 2 



E lab = 2 MeV 




~20 



r, fm 



- Nijmegen—II 

- ISTP, 7h.w 
ISTP, 9hw 



FIG. 22: Small components u p ^{E,r) and Uf^(E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy _Ei ab = 2 MeV. See Fig. ^| for details. 
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FIG. 23: Large components u p ^(E,r) and Uf/f\(E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy _Ei ab = 10 MeV. See Fig. ^| for details. 



FIG. 20: np scattering mixing parameter e in the coupled pf 
waves. See Fig. 1181 for details. 
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FIG. 24: Small components u p /f,(E,r) and Uf^(E,r) of the FIG. 27: Large components u p / p ^(E,r) and Uftf\(E,r) of the 



coupled pf waves np scattering wave function at the labora- 
tory energy i?i ab = 10 MeV. See Fig. ^| for details. 



coupled pf waves np scattering wave function at the labora- 
tory energy _Ei ab = 150 MeV. See Fig. ^|for details. 
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FIG. 25: Large components u p ^(E,r) and u^^(E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy i?i ab = 50 MeV. See Fig.Q^Jfor details. 



FIG. 28: Small components u p ^{E,r) and Uf^(E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy i?i ab = 150 MeV. See Fig. [K]for details. 
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FIG. 26: Small components u p ij\(E,r) and Uj^{E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy i?i ab = 50 MeV. See Fig. ^| for details. 



FIG. 29: Large components u p t p \(E,r) and Uf/f\(E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy £?i ab = 250 MeV. See Fig.[TS]for details. 
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FIG. 30: Small components u p ^{E,r) and u^ p ^(E,r) of the 
coupled pf waves np scattering wave function at the labora- 
tory energy E^ b = 250 MeV. See Fig.[TH]for details. 



TABLE IX: Non-zero matrix elements in hu> units of the Thui 
ISTP matrix in the pf coupled partial waves. 



V pp , matrix elements 



n 

1 
2 
3 



V pp i-i = V pp 

ra.n+l n+l,n 



-0.08320586302194144 
-0.1733874789751888 
-0.1630792532684091 
-0.02514449050450729 



0.06828130087634332 
0.09710466067410284 
0.04737005443342412 



V/J, matrix elements 



V, 



ff 



v, 



ff 



n,n-\-l 



V. 



ff 



n-\- 1 ,n 



-0.01860731179640451 
-0.01230112258503713 
-0.002274165031966646 



0.008146529480927311 
0.002878668408624830 



n 

1 
2 
3 



ypf = yIP 
r n.n—L " n — l,n 



V , matrix elements 

nn' 

fp 



ypf — yfp 

v nn v nn 



-0.02731096515966312 
-0.005320397951202732 
0.009906839670436384 



0.03113837433227350 
0.02654889981453685 
- .007039900977944484 



of Nijmegen-11 at high energies. At the same time, the 
'small' (non-diagonal) ISTP wave function components 
u p rfj(E,r) and Uf< p \{E,r) differ essentially at small dis- 
tances from the Nijmegen-II ones. It is a clear indication 
of a very different nature of the ISTP tensor interaction. 

Now we apply the inverse scattering J-matrix approach 
to the coupled sd partial waves and obtain the 8fku ISTP 
hereafter refered to as Version ISTP. The description of 
the phenomenological data by this potential (and other 
ISTP versions discussed later) is shown in Figs. I^TTOHI 
The np s wave and d wave phase shifts S s and Sd are 
excellently reproduced up to the laboratory energy of 
350 MeV. There is a small discrepancy between the ex- 
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FIG. 31: 3 si np scattering phase shifts 8 S (coupled sd waves). 
Filled circles — experimental data of Ref. Q]; solid line — re- 
alistic meson exchange Nijmegen-II potential Q phase shifts; 
dashed line — Version and Version 1 ISTP phase shifts; 
dotted line — Version 2 ISTP phase shifts. 
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FIG. 32: 3 di np scattering phase shifts Sd (coupled sd waves). 
See Fig. 0U for details. 
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FIG. 33: np scattering mixing parameter e in the coupled sd 
waves. See Fig. OH for details. 
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perimental and the Version ISTP mixing parameter e 
at the laboratory energy of w 25 MeV. However, 
the overall Version ISTP description of experimental 
scattering data (including the mixing parameter e) over 
the full energy interval Ey^ = 0-7- 350 MeV is seen from 
Figs. OH] 031 to be competitive with the Nijmegen-II, one 
of the best realistic meson exchange potentials. 

The Version ISTP is constructed by fitting the ex- 
perimental scattering data, the deuteron ground state en- 
ergy Ed, the s wave asymptotic normalization constant 

s^s and r\ — —-. However, there are other important 
deuteron observables known experimentally such as the 
deuteron rms radius (r 2 ) and the probability of the 
d state. Various deuteron properties obtained with the 
Version ISTP (and other ISTP versions discussed later) 
are compared in Table Ixl with the predictions obtained 
with Nijmegen-II potential and with recent compilations 
of the experimental data |43J, |4J] . It is seen from the ta- 
ble that the Version ISTP overestimates the deuteron 
rms radius and underestimates the d state probability. 

The deuteron wave functions can be calculated by uti- 
lizing the J-matrix formalism at the negative energy 
as is discussed in Ref. 0, . The plots of the deuteron 
wave functions are presented in Fig. 02] It is seen that 
the Version ISTP s wave component is very close to 
that of Nijmegen-II. The Version ISTP d wave compo- 
nent coincides with that of Nijmegen-II at large distances 
since both potentials provide the same s^d value; how- 
ever at the distances less than 5 fm the Version ISTP 
d wave component is suppressed. We note also that the 
Version ISTP scattering wave functions (not shown in 
the figures below) are significantly different from those of 
Nijmegen-II at short distances. 

Our conclusion is that the Version ISTP does not 
seem to be a realistic NN potential. 
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To improve the description of the deuteron properties, 
it appears natural to apply to our Version ISTP a 
phase equivalent transformation that leaves unchanged 
the scattering observables S s , 84, e, the deuteron ground 
state energy Ed and the deuteron asymptotic normaliza- 
tion constants si s and s^d- The phase equivalent transfor- 
mation discussed in Refs. [16l XVA lisj is very convenient 
for our purposes since it is defined in the oscillator ba- 
sis. This transformation gives rise to an ambiguity of 
the potential fit within the inverse scattering J-matrix 
approach, which have been mentioned several times al- 
ready. We now need to discuss this in more detail. 

This phase equivalent transformation is based on the 
unitary transformation 



U= £ E E E K) (n'T'\ 

T=s,d r'=s,d n=0 n'=0 



(66a) 



where the unitary matrix \U] with matrix elements U^, 
should be of the form [H Eld 



[U] = [U ] © [7] = 














(66b) 



and [I] is the infinite unit matrix. The unitary transfor- 
mation 1)66(1 is applied to the infinite Hamiltonian matrix 
[H] in the oscillator basis {|nr)}: 



[H] = [U] [H] [{/+] 



(67) 



The transformed Hamiltonian H is defined through its 
(infinite) matrix [H] with matrix elements HYE, = 

(nT\H\n'r'). That is, the matrix [H] is obtained by 
means of the unitary transformation ((67JI in the original 
basis {\nT)} and not in the transformed basis {\nT}} = 
U{\nT)}. Clearly the spectra of the Hamiltonians H and 
H are identical. If the submatrix [Uq] is small enough, 
the unitary transformation l|67|) leaves unchanged the last 
components (iV r r|A) of the eigenvectors («r|A) obtained 
by solving the algebraic problem 1)40(1 and hence it leaves 
unchanged the functions Qrr' that completely determine 
the fC-matrix, the S'-matrix, the phase shifts S s and 6d, 
the mixing parameter e, the asymptotic normalization 
constants stf s and s^d, etc. 

The potential V entering the Hamiltonian H, phase 
equivalent to the initial potential V entering the Hamil- 
tonian H, can be expressed as 



V = V + AV 7 



(68a) 



r, fm 



where 



FIG. 34: Radial deuteron wave functions. Solid line — realis- 
tic meson exchange Nijmegen-II potential wave functions; 
dot-dash line — Version ISTP wave functions; dashed line 
- Version 1 ISTP wave functions; dotted line — Version 2 
ISTP wave functions. 



AV = H - H. (68b) 

We should improve the tensor component of the NN 
interaction to increase the d state probability in the 
deuteron and reduce the rms radius. Therefore the only 
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TABLE X: Deuteron property predictions obtained with various 8hu> ISTP versions and with Nijmegen-II potential in compar- 
ison with recent compilations |43.l44| . 



Potential 


E d , MeV 


d state rms radius, 
probability, % fm 


St., fm" 1 / 2 


B?d 


Version 


-2.224575 


0.4271 


1.9877 


0.8845 


0.0252 


Version 1 


-2.224575 


5.620 


1.9997 


0.8845 


0.0252 


Version 2 


-2.224575 


5.696 


1.968 


0.8629 


0.0252 


Nijmegen-II 


-2.224575 


5.635 


1.968 


0.8845 


0.0252 


Compilation \43j 


-2.224575(9) 


5.67(11) 


1.9676(10) 


0.8845(8) 


0.0253(2) 


Compilation [44] 


-2.224589 




' 1.9635 
1.9560 
1.950 


0.8781 


0.0272 



non-trivial submatrix [Uq] of the matrix (|66b(l should cou- 
ple the oscillator components \ns) and \n'd) of different 
partial waves. We take the simplest form of the sub- 
matrix [Uq]: a 2 x 2 matrix coupling the |0s) and \0d) 
basis functions. In other words, the non-trivial matrix 
elements U^, constitute a 2 x 2 rotation matrix with a 
single continuous parameter -d: 



[U ] 



U ss U 



TT" 



d 
00 
rjdd 
U Q0 



cos + sin 
— sin d cos i) 



while all the remaining matrix elements 



U, 



it' 



Srr' 



for n > or n > 0. 



(69a) 



(69b) 



Varying the parameter d of the transformation l|67[) - 
(|69|l . we obtain a family of phase equivalent potentials 
and examine which of them provides the better descrip- 
tion of the deuteron properties and np scattering wave 
functions. The best result seems to be the potential ob- 
tained with $ = —14°. This potential is hereafter referred 
to as Version 1 ISTP. 

As a result of the transformation H67 [) - l|69|) . the po- 
tential energy matrix acquires two additional non-zero 
matrix elements \\ff = V$f. These additional matrix 
elements are schematically illustrated by filled circles in 
Fig. 1351 The non-zero matrix elements of the Version 1 
ISTP are given in Table El (in hu = 40 MeV units). 

The deuteron properties obtained with the Version 1 
ISTP are presented in Table Ixl The d state probability is 
improved by the phase equivalent transformation. How- 
ever the phase equivalent transformation produces an in- 
crease of the deuteron rms radius; so this observable be- 
comes even worse than that given by the Version ISTP. 
We found it impossible to obtain an exact description of 
all deuteron properties by means of the phase equivalent 
transformation (|67|l with the simplest matrix i|69|) . 

The deuteron wave functions provided by the Ver- 
sion 1 ISTP are shown in Fig. 03 The Version 1 ISTP 
s wave component is seen to be very close to that of 
the Nijmegen-II. The maximum of the Version 1 ISTP d 
wave component is seen to be shifted to larger distances 





\ SD 







FIG. 35: Structure of the Version 1 and Version 2 ISTP ma- 
trix. The location of non-zero matrix is schematically illus- 
trated by solid lines and filled circles. 



as compared with that of the Nijmegen-II. Of course, the 
shape of the d wave component of the wave function can- 
not be determined experimentally Hence the shape of 
the Version 1 ISTP deuteron wave functions look real- 
istic though these wave functions result in the slightly 
overestimated deuteron rms radius. 

The Version 1 ISTP np scattering wave function com- 
ponents at the laboratory energies E\ a \> = 2, 10, 50, 
150 and 250 MeV are shown in Figs. 1361 1451 in com- 
parison with those of Nijmegen-II potential. As in the 
case of the coupled pf partial waves, the large compo- 
nents u s f s \(E,r) and u d ^(E,r) differ very little from 
the Nijmegen-II ones but the small components are esen- 
tially different at short distances due to the difference of 
the tensor interaction of these two potential models. 

Generally we conclude that the Version 1 ISTP is very 
close to the realistic interaction. The most important 
discrepancy of this interaction is that it overestimates 
the deuteron rms radius by approximately 1.5%. 

We attempted the phase equivalent transforma- 
tion (|67|l with a more complicated matrix [U] than 169fl . 
However, we did not manage to obtain a completely sat- 
isfactory interaction. It is possible to obtain the potential 
providing the required values of the deuteron rms radius 
and of the d state probability by increasing the dimension 
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TABLE XI: Non-zero matrix elements elements in hu units of the Version 1 ISTP matrix in the sd coupled waves. 



n 

1 

2 
3 
4 



n 

1 
2 
3 



V s 



V, 



n 

1 
2 
3 
4 



-0.4576704509056663 
-0.2783240605925930 
-0.01153153008556052 
0.1514476294157645 
0.03632278173824943 



rdd 



v, 



0.008456639591855719 
0.3220439073711110 
0.3084931588662858 

0.06118166034551464 

V-sd T/^ s 
77 . 77. — 1 '77. — 1 .77. 



ti,ti + 



— V s " 

1 — ^ 71+1, 71 



0.2111262515300726 
0.07816883400308394 
-0.05346707187864697 
-0.05592826862748490 



V 71, 71+1 



V 71+1, 71 



-0.08337354364629734 
-0.1788388098603870 
-0.09304409937329216 



V*i — VL/i matrix elements 



v d 

* 77 



T/sd t/^ s 

^71,71+1 — 17 71+1, 71 



-0.06899752955786595 
0.06774418012432615 
0.04913873244946831 

-0.001715094993409672 



-0.4824076895869836 
-0.06136692873982898 
-0.08068524598671112 
-0.02041291263896068 



0.2540123500192352 
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FIG. 36: Large components u 3 ^ a - ) (E,r) and u d ^(E,r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy Ei ah = 2 MeV. See Fig. |M] for details. 



FIG. 37: Small components u a ^ d - ) (E , r) and u d t s \(E,r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy £i ab = 2 MeV. See Fig. |M] for details. 



of the submatrix [Uq] and introducing additional trans- 
formation parameters, but our attempts yielded unreal- 
istic scattering wave functions. 

To improve the sd-ISTP we suggest a slight change to 
the s wave asymptotic normalization constant s# a that 
is used as an input in our inverse scattering approach. 
The stf s value cannot be measured in a direct experi- 
ment. As was mentioned in Ref. |44(> the si a values dis- 
cussed in the literature vary within a broad range from 
0.7592 fm" 1 / 2 to 0.9863 fm" 1 / 2 . Therefore, the modified 
value £/ s = 0.8629 fm -1 / 2 that we use for the construc- 



tion of the improved sci-ISTP, seems to be reasonable. 
We do not change the remaining inputs in our inverse 

scattering approach including rj = — y (and hence we 

modify s^d together with srf s ) to obtain the ISTP of the 
type shown in Fig. 1171 and apply to it the phase equiva- 
lent transformation (|67ll with the parameter d — —14° of 
the matrix (|69|l . This potential is referred to as Version 2 
ISTP. This potential has the structure schematically de- 
picted in Fig. 1351 and its matrix elements are listed in 
Table IXTTl 
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TABLE XII: Non-zero matrix elements elements in hu> units of the Version 2 ISTP matrix in the sd coupled waves. 



n 

1 

2 
3 
4 



n 

1 
2 
3 



n 

1 
2 
3 
4 



V s 



-0.4660631463496376 
-0.2761680294726432 
-0.009473803658917924 
0.1528737342886162 
0.03754792988022171 



rdd 



dd 



0.008667454659207596 
0.3221264718049914 
0.3085166730609980 

0.06120003719298150 

V'sd \/d- s 
n.-Tt. — 1 ' ' ' ii. — 1 . n. 



V, 



1 — v n+l,n 



0.2168839488356998 
0.08090773569137233 
-0.05188144310822707 
-0.05519358984226530 



n,ra+l 



-0.08333937455951757 
-0.1788087936408669 
-0.09301260476557038 



Kf_/ = V„? n matrix elements 



v d 

* n 



\/sd ~[/ds 

v n,n + l — v n + l,n 



-0.06722102540443002 
0.06804449696337271 
0.04940057881591606 

-0.001503998138989182 



-0.4833085003127391 
-0.06047658569273850 
-0.08018710645793066 
-0.02020564623073210 



0.2540038307090694 




FIG. 38: Large components u 3 ^ a - ) (E,r) and u d ^(E,r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy Ei a b = 10 MeV. See Fig. [M] for details. 



FIG. 39: Large components u s t d \(E, r) and u d r s \(E, r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy _Bi ab = 10 MeV. See Fief. I5H for details. 



The deuteron properties are seen from Table H to be 
well described by the Version 2 ISTP. The Version 2 
ISTP scattering wave functions are very close to those 
of the Version 1 ISTP (see Figs. 1351 EESj) . Its deuteron 
wave functions are very close to those of Version 1 ISTP 
(see Fig. 1340 and differ from those of Nijmegen-II in the 
position of the d wave component maximum. 

We suppose that the Version 2 ISTP can be treated as 
a realistic interaction in the coupled sd partial waves. 



IV. APPLICATION OF NN ISTP IN 3 H AND 
4 He CALCULATIONS 

We employ the obtained ISTP in the 3 H and 4 He cal- 
culations within the no-core shell model [III IT2I ] with 
hw = 40 MeV. The same iVTV potentials are used to 
describe the neutron-neutron and neutron-proton inter- 
actions; in the proton-proton case these potentials are 
supplemented by the Coulomb interaction. 

The calculations are performed in the complete NHlj 
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FIG. 40: Large components u 3 ^ a - ) (E,r) and u d ^(E,r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy i?i ab = 50 MeV. See Fig. 01 for details. 



FIG. 43: Small components u a ^ d - ) (E,r) and u d ^(E,r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy E lllb = 150 MeV. See Fig. Olfor details. 
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FIG. 41: Small components u a ( d - ) (E,r) and u d ^{E,r) of the FIG. 44: Large components u s r s \{E, r) and u d r d \(E,r) of the 
coupled sd waves np scattering wave function at the labora- coupled sd waves np scattering wave function at the labora- 
tory energy E Ya . h = 50 MeV. See Fig. 01 for details. tory energy E lllb = 250 MeV. See Fig. 01 for details. 




FIG. 42: Large components u s / a \(E,r) and u d ^(E,r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy E lllb = 150 MeV. See Fig. Olfor details. 



FIG. 45: Small components u a r d \(E, r) and u d t a \(E, r) of the 
coupled sd waves np scattering wave function at the labora- 
tory energy E uh = 250 MeV. See Fig.Olfor details. 
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model spaces with N < 14. We use both 7fiw-ISTP and 
9?kj-ISTP in odd partial waves. The 3 H and 4 He nuclei 
are slightly more bound in the case when we use the 7Hu)- 
ISTP in the odd waves. However, the differences are very 
small: less than 15 keV for 3 H and about 40 keV for 4 Hc. 
The sequence of levels in the 4 He spectrum provided by 
the odd wave 7hu-lSTP and by the odd wave 9foj-ISTP 
is the same but the energies of excited 4 He states are 
shifted down in the case of the odd wave 7Sw-ISTP by 
approximately 100 keV or less. Therefore the deviations 
of the 7fiwTSTP predictions from the experimental odd 
wave scattering data at high enough energies seem to 
produce a negligible effect in the 3 H and 4 He calculations. 
At the same time, 7fiwTSTP has a smaller matrix than 
97kj-ISTP and hence is more convenient in applications. 
Below we present only the results obtained with the 7hu>- 
ISTP in the odd partial waves. 

We have presented various versions of ISTP in the cou- 
pled sd partial waves. The choice of ISTP in other partial 
waves is fixed. Using this fixed set of the non-sd-ISTP in 
combination with the Version M srf-ISTP, we have the set 
of potentials that is refered to as the Version M potential 
model in what follows. 

The 3 H ground state energies E t obtained with the 
Version 1 and the Version 2 potential models in Nhw 
model spaces are presented in Fig.0B]as functions of 1/N. 
It is seen that both potential models provide very similar 
E t values. The convergence of the calculations with N 
appears adequate. The ground state energy E t is seen 
from the figure to be nearly a linear function of 1/N. 
Therefore it is natural to perform a linear extrapolation 
to the infinite Nhw model space, i. e. to the point 1/N = 
0. The linear extrapolation using the two results at the 
highest N-values yields E t rs —8.6 MeV in the Version 1 
potential model and in E t ~ —8.7 MeV in the Version 2 
potential model. 

In Fig. 07| we present the results of the He ground 
state energy E a calculations with the same potential 
models. In the He case we also obtain very similar 
results with the Version 1 and the Version 2 potential 
models. It is interesting that the convergence of the 4 He 
ground state energy is better than that of 3 H. In this 
case the curves connecting the E a values deviate from 
the straight lines. Nevertheless we also perform the lin- 
ear extrapolations of E a (l/N) to infinite N using the E a 
values obtained in 12huj and lAhw calculations and ob- 
tain E a w —26.6 MeV in the Version 1 potential model 
and E a w —27.0 MeV in the Version 2 potential model. 

The quality of the linear extrapolation of -E g . s . (1/iV) 
may be tested in the deuteron calculations. In the 
deuteron case, we know the exact result for the in- 
finite NHlo model space ground state energy Ed = 
—2.244575 MeV obtained by the ^-matrix pole calcu- 
lation with our potentials. The Ed results obtained in 
the Nhw model spaces with N < 14 with the Version 1 
and Version 2 sd-ISTP, are shown in Fig. 08] It is seen 
that Ed(l/N) seems to be a linear function in the in- 
terval 4 < N < 14. The linear extrapolation results 




'o 0,05 0,1 0,15 0,2 0,25 0,3 0,35 0,4 0,45 0,5 

1/N 



FIG. 46: 3 H ground state energy obtained in the Nhco no- 
core shell model calculation vs 1/N. H Version 1 potential 

model; x — Version 2 potential model; dashed line — linear 
extrapolation to the infinite Nhw model space based on the 
last two calculated points ; solid and dash-dot lines are to 
guide the eye. 
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FIG. 47: 4 He ground state energy obtained in the NtkJ no- 
core shell model calculation vs 1/N. See Fig. 1461 for details. 



in Ed ~ —2.5 MeV that differs from the exact energy. 
Therefore the linear extrapolation results can be regarded 
only as a rough estimate of the binding energy. However 
in the 4 He case we achieved a reasonable convergence 
and by the linear extrapolation we increase the binding 
energy by approximately 0.3 MeV only. Therefore our 
estimate of the 4 He binding energy seems to be accurate 
enough. 

The differences in convergence rates for the deuteron, 
3 H and 4 He can be understood from the fact that fkv — 
40 MeV is more optimal for the tighter bound 4 He than 
for the lesser bound systems. 

Our results of the 3 H and 4 He ground state energy cal- 



25 



TABLE XIII: 3 H and 4 He ground state energies (in MeV) obtained in 14ftw no-core shell model calculations and by the linear 
extrapolation to the infinite Nftui model space. 



Potential 
model 


Uhu 


3 H 


extrapolation 


Uhuj 


4 He 


extrapolation 


Version 


-9.091 




9.7 


-33.223 




-33.4 


Version 1 


-7.718 




-8.6 


-26.241 




-26.6 


Version 2 


-7.860 




-8.7 


-26.734 




-27.0 


Nature 




-8.48 






-28.30 






E, MeV 
10 



0,05 0,1 0,15 0,2 0,25 0,3 0,35 0,4 0,45 0,5 

1/N 

FIG. 48: Deuteron ground state energy obtained in the Nhui 
no-core shell model calculation vs 1/N. See Fig. l46l for details. 



culations are summarized in Table IXlTTl We also present 
in the table the results obtained with the less realistic 
Version potential model. Both 3 H and 4 He are essen- 
tially overbound in this potential model. With both Ver- 
sion 1 and Version 2 potential models we obtain a rea- 
sonable description of the 3 H and 4 He bindings. Our 4 He 
results are better than the ones obtained (see Ref. 
with any of the realistic meson exchange interactions 
without allowing for the three-body interactions. In the 
3 H case, we have underbinding in the 14ftu model space 
and a small overbinding obtained by the linear extrap- 
olation. Unfortunately, the difference between the 14Sw 
model space and the linear extrapolation results is rather 
large. Most probably the 3 H ground state energy curve 
in Fig. will flatten out in larger model spaces. This 
will shift the extrapolated ground state energy upwards 
from our current result. Hence the expected ground state 
energy in the N — > oo limit lies between the 14fiw and 
the present linear extrapolation. In other words, our lin- 
ear extrapolation and lAhoj results are expected to be the 
lower and upper boundaries for the exact results, respec- 
tively. An approximately 0.9 MeV difference between 
the lAfkv and the linear extrapolation ground state en- 
ergies in the 3 H case indicates the 0.9 MeV uncertainty 
of our predictions. The 3 H ground state energy obtained 



experiment 



ver.2 





+ 



0+0 ^ 

extrapolation 



10 



-20 



-30 



FIG. 49: 4 He spectrum obtained with Version 2 potential 
model in the no-core shell model in the Uhu (13hui) model 
space for even (odd) parity states. Dashed line shows the 
result of the linear extrapolation of the ground state energy 
to the infinite Nhu model space. Experimental data are taken 
from 451. 



in Faddeev calculations with CD-Bonn NN potential is 
—8.012 MeV (see |f|). All the remaining modern realistic 
meson exchange potentials predict the 3 H binding energy 
to be less than 7.4 MeV 0- Therefore our 3 H binding en- 
ergy predictions are not worse than those obtained with 
the realistic meson exchange potentials without allowing 
for the three-body forces while our 4 He binding energy 
predictions are better. 

In Fig. 02| we present the spectrum of the lowest ex- 
cited He states of each J 7r '. The description of the ex- 
cited states energies is reasonable though further from 
experiment than the ground state. On the other hand, 
we expect the excited states to be less converged and to 
drop more in larger model spaces. Of course, a full dis- 
cussion of the states above breakup must await proper 
extensions of the theory to the scattering domain. 



V. CONCLUDING REMARKS 

We obtained nucleon-nucleon ISTP potentials by 
means of the J-matrix version of the inverse scattering 
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approach. The potentials accurately describe the scatter- 
ing data. They are in the form of 8ftu;-truncated matrices 
in the oscillator basis with hui = 40 MeV. The potential 
matrices are tridiagonal in the uncoupled partial waves. 
In the coupled partial waves, the potential matrices have 
two additional quasi-diagonals in each of the submatri- 
ces responsible for the channel coupling. The sd-ISTP of 
this type (Version 0) underestimates the deuteron d state 
probability and overestimates the deuteron rms radius. 
We designed two other sd-ISTP with two additional ma- 
trix elements providing the correct description of the d 
state probability, one of them (Version 1) overestimates 
the rms radius by approximately 1.5% while the other 
one (Version 2) provides the correct description of the 
deuteron rms radius. All other deuteron observables are 
reproduced by all sd-ISTP versions. 

The ISTP potentials are used in the 3 H and 4 He no- 
core shell model calculations. Both Version 1 and Ver- 
sion 2 ISTP potential models provide very good predic- 
tions for the 3 H and 4 He binding energies and a rea- 
sonable 4 He spectrum. With the less realistic Version 
potential model, we obtain overbound 3 H and 4 He nu- 
clei. We note that there were other attempts to design 
the NN interaction providing the description of the tri- 
ton binding energy together with the NN scattering data 
and the deuteron properties Our interactions are 

much simplier and can be directly used in the shell model 
calculations of heavier nuclei. 

Generally our approach is aimed at shell model appli- 
cations in heavier nuclei. However our potentials are sim- 
ple enough and can be used directly in other microscopic 
approaches, e. g. in Faddeev calculations. We hope that 
our interactions minimize the need for three-body forces. 
It is known ^(j that the three-body force effect can be re- 
produced in a three-body system by the phase equivalent 
transformation of the two-body interaction. This phase 
equivalent transformation can also spoil the description 
of the deuteron observables, in particular, the deuteron 
rms radius can be arbitrary changed by phase equivalent 
transformations • We expect that there exist transfor- 
mations minimizing the need for three-body force effects, 



that do not significantly change the nucleon-nucleon in- 
teraction. That is, the deuteron properties, the deuteron 
and scattering wave functions of the transformed NN po- 
tential may remain very close to the ones developed here 
while achieving improved descriptions of other nuclei. In 
this context, it is worth noting that our approach does 
not assume cither a particular operator structure to the 
interaction or locality. 

From this point of view, the Version 2 ISTP accurately 
describing the deuteron properties and providing good 
predictions for the 3 H and 4 He bindings, can be regarded 
as such an interaction effectively accounting for effects 
that might otherwise be attributed to three-body forces. 
Clearly, additional efforts may provide superior NN in- 
teractions with less dependence on three-body forces for 
precision agreement with experiment. 

Finally, we suggested a new approach to the construc- 
tion of the high-quality NN interaction and examined 
the obtained ISTP NN interaction in three and four nu- 
cleon systems by means of the no-core shell model. The 
3 H and 4 He binding energies are surprisingly well de- 
scribed. Obviousely it will be very interesting to extend 
these studies on heavier nuclei, to investigate in detail 
not only their binding but the spectra of excited states 
as well. It is also important to investigate more carefully 
the ISTP description of the two-nucleon system since, for 
example, we have deferred the discussion of the deuteron 
quadrupole moment Q. We just mention here that the 
Version 2 ISTP prediction of Q = 0.317 fm 2 is not so far 
from the experimental value of 0.2875 ± 20 fm 2 0. The 
phase equivalent transformations discussed above make 
it possible to improve the Q predictions and to examine 
the effect of such improvement in light nuclear systems. 
We plan to address this problem in future publications. 

This work was supported in part by the State Pro- 
gram "Russian Universities" , by the Russian Foundation 
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